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O ! Abstract 

o 

i/-^ ■ The renormalization of the Schrodinger equation with regular One Boson Ex- 

change and singular chiral potentials including One and Two-Pion exchanges 
is analyzed within the context of NN scattering. 

^ : 1 Introduction 

One traditional view of NN force has been through One Boson Exchange 
(OBE) Models [1,2]. Recent developments have shown how chiral symmetry 
K*" \ may provide NN forces of practical interest in nuclear physics [3-5] . Remark- 

ably, chiral expansions, based on assuming a large scale suppression on the 



O 



■ parameters Airf-n ~ Mn ~ IGeV necessarily involve singular potentials at 

^ . short distances, i.e. r^|V"(r)| — »• cxo for r — > 0. If we take the limit r -C l/rnvr 

^ ! (or equivalently large momenta) pion mass effects are irrelevant and hence 

\ at some fixed order of the expansion one has 

; ^ (47r/,)2«M™ r2«+'^ 

H \ 

I (the only exception is the singlet channel-OPE case which behaves as ~ 

ml/ f^r, see below). The dimensional argument is reproduced by loop cal- 
culations in the so called Weinberg dimensional power counting [6,7]. Thus, 
much of our understanding on the physics deduced from chiral potentials 
might be related to a proper interpretation of these highly singular poten- 
tials. Renormalization is the most natural tool provided 1) we expect the 
potential is realistic at long distances and 2) we want short distance details 
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not to be essential in the description. This is precisely the situation we 
face most often in nuclear physics. Knowledge on the attractive or repulsive 
character of the singularity turns out to be crucial to successfully achieve this 
program and ultimately depends on the particular scheme or power counting 
used to compute the potential. We illustrate our points for the simpler OBE 
potential in the ^5*0 channel and then review some results for chiral OPE and 
TPE potentials for all partial waves and the deuteron bound state. 



2 Renormalization of OBE potentials 

The singularity of chiral potentials raises suspicions and, quite often, much 
confusion. However, if properly interpreted and handled they do not dif- 
fer much from the standard well-behaved regular potentials one usually en- 
counters in nuclear physics. Actually, we digress here that renormalization 
may provide useful insights even if the potential is not singular at the origin 
{r'^V{r) !). For definiteness, let us analyze as an illustrative example the 
phenomenologically successful ^5*0 OBE potential [1,2] (we take rup = m^) 

Y(j^) = yTrNN"''n ^ UaNN ^ , iiuNN ^ , /cys 

^ ' ~ IQi^Ml r An r Att r ^ ' 

where for simplicity we neglect nucleon mass effects and a tiny rj contribution. 
We take = 138MeV, Mn = 939MeV, = 783MeV and g^NN = 13.1 
which seem firmly established. Actually, Eq. looks like a long distance 
expansion of the potential. NN scattering in the elastic region below pion 
production threshold involves CM momenta p < Pmax = 400MeV. Given 
the fact that l/m^ = 0.25fm ^ 1/pmax = 0.5fm we expect heavier mesons 
to be irrelevant, and u itself to be marginally important. In the traditional 
approach, however, this is not so [1,2]. Actually, the problem is essentially 
handled by solving the Schrodinger equation (S-wave) 

- u'lir) + MVir)up{r) = p\ir) (3) 

with the regular solution at the origin, Up{0) = 0. This boundary condition 
implicitly assumes taking also the potential all the way down to the origin. 
The asymptotic condition for r ^ l/wijr is taken to be 

sin(pr + (5o(p)) 
sm6o{p) 

where 6o{p) is the phase-shift. For the potential in Eq. ([21) the phase shift is 
an analytic function of p with the closest branch cut located at p = ±im.„^/2, 
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so that one can undertake an effective range expansion, 

p cot 5o{p) = — ^ + \rQp^ + V2p'^ + . . . (5) 
tto ^ 

within a radius of convergence \p\ < 171^^/2. A similar expansion for the wave 
function Up{r) = Uo{r) +p^-U2(r) . . . means solving the set of equations 

- <(r) + M\/(r)no(r) = 0, (6) 
Uo(r) 1-r/ao, 

-u'^{r) +U{r)u2{r) = Uo{r) , (7) 
U2{r) (r^ - 3aor^ + Saoror) / (6ao) , 

where, again, the regular solutions, uo{0) = U2{0) = are taken. With this 
normalization the effective range tq is computed from the standard formula 

/■oo 

ro = 2 / dr [(1 - r/aof - Mo(r)'] . (8) 

In the usual approach [1,2] everything is obtained from the potential assumed 
to be valid for < r < oo. In practice, strong form factors are included 
mimicking the finite nucleon size and reducing the short distance repulsion 
of the potential, but the regular boundary condition is always kept. [f| As it is 
well known the ^5*0 scattering length is unnaturally large = — 23.74(2)fm, 
while ro = 2.77(4)fm. Let us assume we have fitted the potential, Eq. ([2]), 
to reproduce ao- Under these circumstances a tiny change in the potential 
V ^ V + AV has a dramatic effect on ao, since one obtains 

POO 

Aao = alMN AV{r)uo{r)^dr . (9) 
Jo 

As a result, potential parameters must be fine tuned. In particular, the 
resulting cu-repulsive contribution is well determined with an unnaturally 
large coupling, g^NN ~ 16. [1,2]. In our case, with no form factors nor 
relativistic corrections, a fit to Ref. [8] yields g^jNN = 12.876(2), gaNN = 
12.965(2) and = 554.0(4)MeV with xVDOF = 0.26. Note the small 
uncertainties, as expected from our discussion and Eq. (Q. As mentioned 
above 1/m^ = 0.25fm -C 1/pmax = 0.5fm so u should not be crucial at least 
for CM momenta p <^ Pmax- Thus, despite the undeniable success in fitting 
the data this sensitivity to short distances looks counterintuitive. 



■^Calculations solving the equivalent Lippmann-Schwinger equation in momentum space 
for regular potentials correspond always to choose the regular solution for the Schrodinger 
equation in coordinate space. 
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The renormalization viewpoint refuses to access physically the very short 
distance region, but encodes it through low energy parameters described by 
the effective range expansion, Eq. as renormalization conditions (RC's). 
In the case of only one RC where ao is fixed one proceeds as follows [9, 10]: 

• For a given ao integrate in the zero energy wave function uo{r), Eq. ([6]), 
down to the cut-off radius r^. This is the RC. 

• Implement self-adjointness through the boundary condition 

Up{rc)uo{rc) - Mo(rc)up(rc) = , (10) 

• Integrate out the finite energy wave function Up{r), from Eq. ([3]) and 
determine the phase shift So{p) from Eq.(jl]). 

• Remove the cut-off — to strive for model independence. 

This allows to compute 6o (and hence ro, W2 ) from V{r) and ao as independent 
information. Note that this is equivalent to consider, in addition to the 
regular solution, the irregular one □. A beautiful result is the universal low 
energy theorem which highlights this de-correlation between the potential 
and the scattering length [10] 

ro = 2 / c/r(l-M^o) / c/r(r - Mo.oMo.i) + ^ / dr {r"^ - ul (11) 

Jo ' ao Jo ' ag Jo 

based on the superposition principle of boundary conditions, i.e. writing 
Uo{r) = Uofi{r) — Mo,i('")/«o with Uo,n{^) f"' and using Eq.([H]). A fit of 
the potential ^ with q^jnn = to the effective range yields (Fig. [T]) a 
strong correlation between mg- and QaNN- Over- imposing this correlation to 
ro = 2.670(4)fm, a fit to Ref. [8] yields = 493(12)MeV, ^^^iv = 8.8(2), 
gujNN = 0(5) with x^/DOF = 0.24 (Fig. [T]). Note the larger uncertainties, 
although correlations allow Qujnn ~ 9 and m„ ~ 520MeV within Ax^ = 1- 
Contrary to common wisdom, but according to our naive expectations, no 
strong short range repulsion is essential. The moral is that building ao from 
the potential is equivalent to absolute knowledge at short distances and in 
the ^5*0 channel a strong fine tuning is at work. Of course, a more systematic 
analysis should be pursued in all partial waves and relativistic corrections 
might be included as well, but this example illustrates our point that the 
renormalization viewpoint may tell us to what extent short distance physics 
may be less well determined than the traditional approach assumes. This 
opens up a new perspective to the phenomenology of OBE potentials in cases 
where the strong a;-repulsion has proven to be crucial at low energies [12]. 

^In momentum space this can be shown to be equivalent to introduce one counterterm 
in the cut-ofF Lippmann-Schwinger equation, see Ref. [11] for a detailed discussion. 
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hIq - correlalion 



Re normalized neutron-proton phase shifts 



o^^- 23.1 A 
r,, (ejcp) - 2.77 (5) 
500 - r|,{Niimegen)^2.67 




(dimensionless) 
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- gaNN = 8.8(2) 




8o)NN = (5) 




K^IDOP^ 0.26 





200 250 
k [MeV] 



Figure 1: Results in the ^5*0 channel for the renormalized OBE potential. 
Left: Effective range correlation between g^NN and mo- for Qujnn = 0. Right: 
Renormalized phase shift (in degrees) as a function of the CM momentum 
(in MeV) in the OBE {it + a + u) model. The data are an average of [8]. 



3 Renormalization of chiral potentials 

The generalization of the above method to the singular chiral potentials [6,7] 
has been implemented in [10] with promising results for One- and Two Pion 
Exchange (OPE and TPE). We illustrate again the case of pn scattering in the 
^5'o-channel. For the simplest situation with one RC, corresponding to fix the 
scattering length as an independent parameter, the method outlined above 
may be directly applied to singular potentials provided they are attractive, 
i.e. V{r) —Cn/r^ with n > 20- The result for zero energy wave functions 
as well as the effective range can be seen at Fig. [21 NNLO corresponds to 
the TPE potential of Ref. [6]. As we see the Nijmegen result ro = 2.67fm 
is almost saturated by the TPE potential yielding tq = 2.87fm already at 
Tc ~ 0.5fm. Calculations with TPE to N3L0 with one RC show convergence 
but no improvement [11] without or with A explicit degrees of freedom. Thus, 
some physics is missing, perhaps Svr effects. If, in addition to ao, we want to 
fix To = 2.67fm [8] as a RC we must solve Eqs. (Ej) and (171). The matching 
condition at the boundary r = becomes energy dependent [13] 

K^^c) ^ M^(re)+pV2(r,) + ... 
Up{rc) Uo{rc) + p^U2{rc) + . . . ' 

The generalization to arbitrary order is straightforward. For RC's we have 
'^pi'^) — ^n=o^'^"^2n('") and using the natural extension of the matching rela- 
tion in Eq. ( lT2l) as well as the superposition principle of boundary conditions 

^If the potential was singular and repulsive one cannot fix any low energy parameters; 
doing so yields non-converging phase shifts. 
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one can show the following formula 

pcot5o(p) = gr^^47^' ^13) 

where the coefficients be related to the effective range parameters 

= 1, «! = — l/ao, (32 = rQ,a3 = V2 etc. and An{p,rc) and Bn{p,rc) are 
functions which are finite in the limit Tc — > and depend solely on the 
potential. In Eq. (fT3l) the dependence on the low energy parameters used as 
input is displayed explicitly and can be completely separated from the long 
range potential [13]. The coupled channel case can be analyzed in terms of 
eigenpotentials although the result is cumbersome. In Fig. [S] we show the 
phase shitf for the ^5*0 channel when the potential is considered at LO, NLO 
and NNLO and either one RC (fixing oq) or two RC's (fixing ao and ro) are 
considered. LO+IC, NL0+2C and NNL0+2C fix the same number or RC's 
as LO, NLO and NNLO of the Weinberg counting respectively. As we see, 
our NNL0+2C does not improve over NL0+2C. 

It is worth mentioning that the innocent-looking energy dependent match- 
ing condition, Eq. (fT2l) . is quite unique since this is the only representation 
guaranteeing finiteness of results for singular potentials [13]. Polynomial ex- 
pansions in such as suggested e.g. in Ref. [7] do not work for Vc ^ 0. A 
virtue of the coordinate over momentum space is that these results can be de- 
duced analytically. For instance, the equivalent representation of Eq. ffT^ in 
momentum space may likely exist, but is so far unknown. Actually, the usual 
polynomial representation of short distance interactions in momentum space 
Vs{k\ k) = Co + C2(A;2 + k'^) + ... of standard NLO and NNLO Weinberg 
counting is renormalizable only when C2 — for A — oo [11]. 



4 Renormalization of the Deuteron 



In the ^Si—^Di channel, the relative proton-neutron state for negative energy 
is described by the coupled equations 



-|^ + M^K(r) 
MNVsdir) 



-7 



(14) 



Here 7 = \/MmB, with B = 2.24MeV is the deuteron binding energy and 
u{r) and w{r) are S- and D-wave reduced wave functions respectively. At 
long distances they satisfy. 



Ase 



-■yr 



1 

7r (7^")^ 



(15) 
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Figure 2: (Left panel) Zero-energy, ^5*0 linearly independent wave functions 
at NNLO; Mi — 1 and Ur ^ r for r oo. (Right Panel) Effective range tq 
as a function to the cut-off for the same channel and different orders; using 
ro(re) = 2 (/~(1 - r/«o)' dr - /~ c/r) , with «o = -23.74fm [8]. 



where rj is the asymptotic D/S ratio parameter and As is the asymptotic 
normalization factor, which is such that the deuteron wave functions are 
normalized to unity. The OPE ^Si Di potential is given by Mi^Vg = Uc 
, Mj^Vsd = 2V2UT,MNVd = Uc- 2Ut where for r > > we have 



^ ~ 16^7/2 ' I^T- Uc[l + — + ) . (lb) 

The tensor force generates a singularity at the origin in coupled chan- 
nel space. This behavior of the potential is strong enough to overcome the 
centrifugal barrier at short distances, thus modifying the usual threshold be- 
havior of the wave functions. The interesting aspect of this potential is that 
after diagonalization it has one positive (repulsive) and one (negative) attrac- 
tive eigenvalue. The proper normalization of the wave functions in the limit 
rc —>■ implies that one can only fix one free parameter, e.g. the deuteron 
binding energy [10]. Other properties may be predicted, for instance one 
gets r]oPE = 0.0263 (exp. 0.0256(4)). The TPE chiral potentials of Ref. [6] 
have also been renormalized [10], yielding a rather satisfactory picture of the 
deuteron. The results described here have been reproduced in momentum 
space [14] . The required cut-off in momentum space is larger than a naive es- 
timate ~ 1/A because the regularization influences both the counterterms 
as well as the potential. Deuteron form factors, probing some off-shellness of 
the potential, have been computed describing surprisingly well the data up 
to momenta q ~ SOOMeV when LO currents are considered [§. 



®See talk of D. R. Phillips in this conference 
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Figure 3: Renormalized ^5*0 phase shifts (in degrees) for chiral LO,NLO and 
NNLO potentials fixing ao = — 23.74fm (Left panel) or oq = — 23.74fm and 
tq = 2.77fm(Right panel) as input parameters. The data are from [8]. 



5 Power counting and renormalization 

The question on how a sensible hierarchy for NN interactions should be orga- 
nized remains so far open, because it is not obvious if one should renormalize 
or not and how [10,15-17]. However, for a given long distance potential, we 
know whether and, in positive case, how this can be made compatible with 
the desired short distance insensitivity [10]. Not all chiral interactions fit into 
this scheme, and thus it is sometimes preferred to keep finite cut-offs despite 
results being often strongly dependent on the choice at scales Vc ~ 0.5 — Ifm 
similar to the ones we want to probe in NN scattering [10,15]. Renormaliz- 
ability of chiral potentials developing a singularity such as Eq. ([T]) requires 
that one must choose the regular solution in which case the wave function 
behaves as Up{r) ~ (r47r/,r)^^ and thus increasing insensitivity is guaran- 
teed as the power of the singularity increases. Converging renormalized TPE 
calculations show insensitivity for reasonable scales of ~ 0.5fm [10]. 

The Weinberg counting based in a heavy baryon approach at LO [5] for 
^5*0 and ^Si—^Di states turns out to be renormalizable. There is at present no 
necessity argument why this ought to be so, for the simple reason that power 
counting does not anticipate the sign of the interaction at short distances. 
When one goes to NLO the short distance 1/r^ singular repulsive character 
of the potential makes the deuteron unbound [10]. Finally, NNLO potentials 
diverge as — 1/r^ and are, again, compatible with Weinberg counting in the 
deuteron [10]. More failures have been reported in Refs. [11,15]. Relativistic 
potentials subjected to different power counting have been renormalized in 
Ref. [18] yielding much less counterterms due to their different short distance 
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singularities and slightly better overall description, although the ^5*0 
phase is not improved as compared to the heavy baryon formulation. These 
complications in the more fundamental chiral potentials contrast with the 
simplicity of the a + tt OBE renormalized results (see Figs. [HandlS]). 

In the present state of affairs a clue might come from a remarkable anal- 
ogy between the NN interaction in the chiral quark model and the Van der 
Waals molecular interactions in the Born-Oppenheimer approximation [10]. 
For non-relativistic constituent quarks the direct NN interaction is provided 
by the convoluted OPE quark-quark potential. Second order perturbation 
theory in OPE among quarks generates TPE between nucleons yielding 



When HH' = NA and HH' = AA this resembles Ref. [19] which for 
2fm < r < 3fm behaves as a exchange with = 550MeV and q^nn = 9.4. 
Moreover, the second order perturbative character suggests that the poten- 
tial becomes singular ~ 1/r^ and attractive, necessarily being renormalizable 
with an arbitrary number of counterterms through energy dependent bound- 
ary conditions [13]. Clearly, the renormalization of such a scheme where the 
A^A splitting is treated as a small scale deserves further investigation [20]. 

6 Conclusion 

Renormalization is the mathematical implementation of the appealing phys- 
ical requirement of short distance insensitivity and hence a convenient tool 
to search for model independent results. In a non-perturbative setup such as 
the NN problem, renormalization imposes rather tight constraints on the in- 
terplay between the unknown short distance physics and the perturbatively 
computable long distance interactions. This viewpoint provides useful in- 
sights and it is within such a framework that we envisage a systematic and 
model independent description of the NN force based on chiral interactions. 
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Vnn = {NN\Vope\NN) + J2 



HH'^NN 




+ ... 



(17) 
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